The deterministic nature of conservative post-Newtonian accurate dynamics of 
compact binaries with leading order spin-orbit interaction 
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We formally show that the conservative second post-Newtonian (PN) accurate dynamics of spin- 
ning compact binaries moving in eccentric orbits, when spin effects are restricted to the leading 
order spin-orbit interaction cannot be chaotic for the following two distinct cases: (i) the binary 
consists of compact objects of arbitrary mass, where only one of them is spinning with an arbitrary 
spin and (ii) the binary consists of equal mass compact objects, having two arbitrary spins. We 
rest our arguments on the recent determination of PN accurate Keplerian-type parametric solutions 
to the above cases, indicating that the PN accurate dynamics is integrable in these two situations. 
We compare predictions of our case (i) with those from a numerical investigation of an equivalent 
scenario that observed chaos in the associated dynamics. We also present possible reasons for the 
discrepancies. 

PACS numbers: 04.25.Nx, 95.10.Eg, 97.60.Jd, 97.60.Lf 



I. INTRODUCTION 

Inspiralling binaries consisting of spinning compact ob- 
jects of arbitrary mass ratio are astrophysically interest- 
ing systems. These binaries are the most promising tar- 
gets for the first direct detection of gravitational radia- 
tion by ground-based laser interferometric gravitational- 
wave detectors [1]. The radio- wave detections of spin- 
ning compact binaries are heavily sought after as they 
provide unique laboratories for relativistic gravity and 
plasma physics [2]. The orbital dynamics of such com- 
pact binaries are well described by the post-Newtonian 
(PN) approximation to general relativity. The PN ap- 
proximation to general relativity allows one to express 
the equations of motion for a compact binary as cor- 
rections to Newtonian equations of motion in powers of 
(v/c) 2 ~ GM/(c 2 R), where v, M, and R are the char- 
acteristic orbital velocity, the total mass and the typi- 
cal orbital separation of the binary. To extract lots of 
astrophysical informations from these binaries, both us- 
ing gravitational and radio wave observations, it is un- 
derstood that the compact binary dynamics should be 
known at least to second post-Newtonian (2PN) order 
[3, 4]. The 2PN order gives (v/c) 4 corrections to New- 
tonian equations of motion and the orbital dynamics is 
still conservative at this order. In the case of spinning 
compact binaries, the dynamics is determined, not only 
by the orbital equations of motion for these objects, but 
also by the precessional equations for the orbital plane 
and the spin vectors themselves [5]. The dynamics of 
these systems are required to be deterministic, at least 
during the observational window, for making any mean- 
ingful astrophysical measurements. 

Few years ago, it was argued that spin-orbit and spin- 
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spin contributions to the conservative PN dynamics will 
force the orbits of compact binaries to become chaotic 
[6]. The "chaos" in the conservative dynamics of spin- 
ning compact binaries were identified with the aid of 
phase-space trajectories that form fractal basin bound- 
aries. These trajectories were computed by numerically 
integrating 2PN accurate equations, describing the com- 
pact binary dynamics in harmonic gauge, available in 
Ref. [7]. It was reasoned that compact binaries, having 
initial conditions near the fractal basin boundaries will 
have unpredictable evolution. Therefore, it was argued 
that spinning compact binaries, evolving under PN ac- 
curate orbital dynamics, will have unpredictable gravita- 
tional waveforms during the inspiral. This will make the 
deployment of "matched filtering" and hence the detec- 
tion of these binaries, via gravitational radiation, an im- 
possible task. The above claim was questioned in Ref. [8]. 
These authors computed the divergence of nearby phase 
space trajectories of spinning compact binaries, with the 
aid of Lyapunov exponents, evolving under the same PN 
dynamics, as investigated in Ref. [6] and found no ev- 
idence for any chaotic behavior. However, the above 
analysis was strongly criticized in Ref. [9] , based on argu- 
ments that in Ref. [8] the "proper" Lyapunov exponents 
were not computed. Further, the subtleties involved in 
the use of Lyapunov exponents in general relativity were 
emphasized in Ref. [10]. A recent detailed investigation 
of the conservative PN dynamics of spinning compact bi- 
naries, to substantiate the results of Ref. [6], using again 
the method of fractal basin boundaries, also observed the 
existence of "chaos" in these systems [11]. Moreover, it 
was even shown that 2PN accurate conservative dynam- 
ics of spinning compact binaries can be chaotic even if 
only one of the compact objects spins [see Figs. 4 and 5 
in Ref. [11]]. 

In this paper, we show that the conservative 2PN ac- 
curate dynamics of spinning compact binaries cannot be 
chaotic, when only one of the compact objects spins. We 
rest our arguments on the fact that the above dynam- 
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ics is integrable. In fact, it can be shown that even the 
conservative 3PN accurate dynamics of spinning compact 
binaries, when spin effects are restricted to the leading 
order spin-orbit interaction, is integrable in two distinct 
cases: (i) the binary consists of compact objects of ar- 
bitrary mass, where only one of them is spinning with 
an arbitrary spin and (ii) the binary consists of equal 
mass compact objects, having two arbitrary spins. The 
integrability is evident from the recent determination of 
PN accurate Keplerian-type parametric solution to the 
above two cases [12, 13]. Our 2PN accurate analytic so- 
lutions, especially for case (i), imply that the associated 
dynamics cannot be chaotic. However, numerical inves- 
tigations did report chaos in 2PN accurate dynamics of 
spinning compact binaries when only one of the compact 
objects spins, i.e., our case (i) [11]. We present possible 
arguments for the above contradiction. 

In what follows, we briefly describe the recently de- 
termined PN accurate parametric solutions and list their 
salient features. We present our arguments why the as- 
sociated dynamics cannot be chaotic. We also clarify 
few erroneous statements made in Ref. [11]. Finally, we 
suggest a way to analyze PN accurate dynamics of spin- 
ning compact binaries, semi-analytically, but in a PN 
accurate manner. This should be helpful to determine if 
inspiralling spinning compact binaries, whose dynamics 
is correctly given by PN approximation to general rel- 
ativity, can exhibit chaotic behavior, for the cases not 
considered here, in a realistic manner. 



II. THE RELEVANT BINARY DYNAMICS AND 
ITS DETERMINISTIC NATURE 

Let us first consider the 2PN accurate conservative dy- 
namics of spinning compact binaries, when the spin ef- 
fects are restricted to the leading order spin-orbit inter- 
action. The dynamics is fully specified by a PN accu- 
rate (reduced) Hamiltonian H , which may be symboli- 
cally written as 

H(r,p,Si,S 2 ) = H N (r,p) + H 1PN (r,p) + H 2PN (r,p) 
+ H S o(r,p,SuS 2 ), (1) 

where Hn, Hipn, and i?2PN are, respectively, the Newto- 
nian, first, and second PN contributions to the conserva- 
tive dynamics of compact binaries, when the spin effects 
are neglected. The leading order spin-orbit coupling to 
the binary dynamics is given by Hgo . In the above equa- 
tion, r = IZj (GAI), r = \r\, and p = V j [i, where TZ and 
V are the relative separation vector and its conjugate 
momentum vector, respectively. The familiar symbols 
M and \x have the usual meaning, namely, the total mass 
and the reduced mass. The explicit expressions for PN 
corrections, associated with the motion of nonspinning 
compact binaries, were obtained in Refs. [14, 15]. The 
spin-orbit contributions, available in Refs. [13, 15, 16], 



read 

H S o(r,p,S l ,S 2 ) = -^(rxp)-S eS , (2) 

where S e s gives the effective spin defined by 

S e « = $iSi + 5 2 S 2 ■ (3) 

In the above equation, Si = ^ + f (l — y 7 ! — 4m an d ^2 = 
2 + | (l + \/l — 4?y) , where r\ is the finite mass ratio r\ = 
[j,/M. The reduced spin vectors Si and S 2 are related 
to the individual spins Si and S 2 by S\ = Si/(fiGM) 
and S 2 = S 2 j (pGM), respectively. We recall that TZ, 
"P, Si, and S 2 are canonical variables, such that the 
orbital variables commute with the spin variables, e.g., 
see Refs. [5, 16]. 

Recently, as mentioned earlier, the 3PN accurate 
Keplerian-type parametric solution to the dynamics of 
spinning compact binaries, when spin effects are re- 
stricted to the leading order spin-orbit interaction, for 
the following two distinct cases, namely, (i) the binary 
consists of compact objects of arbitrary mass, where only 
one of them is spinning with an arbitrary spin and (ii) 
the binary consists of equal mass compact objects, hav- 
ing two arbitrary spins, were obtained in Refs. [12, 13]. 
We display below the 2PN accurate version of the above 
parametric solution. It turned out that PN accurate 
dynamics, we are interested in, allows Keplerian-type 
parametrization in an orbital orthonormal triad, k). 
In this triad, i defines the line of nodes, associated with 
the intersection of the orbital plane with the invariable 
plane (ex,ey), which is the plane perpendicular to the 
total (reduced) angular momentum J = Jez, the only 
constant vector in our dynamics. The unit vector k is 
always perpendicular to the orbital plane and defined by 
k = L/L, where L = r x p is the (reduced) orbital an- 
gular momentum vector. The dynamical vectors of the 
problem are given by 

r(t) = r(t) cos ip(t)i(t) + r(t) sin <p(t) j(t) , (4) 
L{t) = Lk(t) , (5) 
S(t) = Je z - Lk(t) , (6) 

where S is the (reduced) total spin, which is related to 
the effective spin by S e s — XsoS for the considered cases. 
The mass dependent coupling constant x so equals 5i or S 2 
for the single-spin case and is given by Si = 5 2 = 7/8 for 
the equal-mass case. The magnitude of J is denoted by 
J and is given by J = (L 2 + S 2 + 2LScosa) 1/2 . Note 
that the angle a between L and S can be chosen quite 
arbitrarily. The time dependent basic vectors k) are 
explicitly given by 

i — cos Tex + sinTey , (7a) 
j = — cos 9 sin Tex + cos 9 cos Tey + sin 9 ez , (7b) 
k = sin 9 sin Tex — sin 9 cos Tey + cos 9 ez ■ (7c) 

The angle 9 gives the constant precessional angle be- 
tween L and J [see Fig. 1 in Ref. [13]]. The time evo- 
lution for r, ip, and T is given in a parametric and PN 
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accurate way, which reads 

r = a r (1 — e r cos it) , 

; _ / , . \ • . 9it , N . fit . 

I = n (t — to ) = u — et sin u H — — (v — u) H — -r sin v 

c 4 c 4 



(8a) 

'J, 
(8b) 



<y9 — </?o = (1 + k)v + -4r sin 2w + sin 3v . 

c 4 c 4 



T - T r 



XsoJ 



c 2 L 3 

where v = 2 arctan 



(v + esinii) , 
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tan — 
2 



(8c) 
(8d) 

(8e) 



The PN accurate orbital elements ei , k, and 

e^,, and the PN order orbital functions 1744, fa, f<i^, and 
g4 V , expressible in terms of E, L, S, 77, and a, are ob- 
tainable from Ref. [13]. We recall that the above para- 
metric solution is usually referred to as the generalized 
quasi-Keplerian parametrization for the PN accurate or- 
bital dynamics of compact binaries and extends the sem- 
inal works done by Damour, Deruelle, Schafer, and Wex 
[15, 17-19]. The quasi-Keplerian parametric solution is 
heavily employed to construct the "timing formula" for 
relativistic binary pulsars, which in turn is used to ex- 
tract astrophysical information from these systems and 
to test general relativity in strong field regimes [20, 21]. 
The recent determination of "ready to use" search tem- 
plates for nonspinning compact binaries moving in inspi- 
ralling eccentric orbits, given in Ref. [22], also required 
the above mentioned quasi-Keplerian parametrization. 

The existence of such a PN accurate parametric solu- 
tion and hence the integrability of the spinning binary 
dynamics should not be surprising if one carefully counts 
the degrees of freedom and the existing conserved quan- 
tities in the above mentioned two distinct cases, where 
our parametrization is valid. In our Hamiltonian sys- 
tem, for the two cases considered, the dimension of the 
associated phase space is 8 [the nonspinning part of the 
Hamiltonian contributes 6 dimensions and the spin con- 
tributes 2 dimensions. This indicates that the number of 
independent degrees of freedom, for the two cases con- 
sidered, is 4 and there are 4 independent constants of 
motion, namely, E, L, S, and L ■ S. The deterministic 
nature of the Hamiltonian, given by Eq. (1) and treated 
as a self-contained dynamical system, naturally follows 
from above arguments. 

One may argue, that the equations of motion for spin- 
ning compact binaries presented in Ref. [7] and employed 
in Refs. [6, 11] are in harmonic coordinates and employ 
a different spin supplementary condition (SSC) than the 
one used in Ref. [13]. Therefore it is not desirable to com- 
pare results of Ref. [11] with results of Ref. [13], especially 
about the deterministic nature of the underlying dynam- 
ics. Following Ref. [12], it is trivial to obtain a harmonic 
coordinate version of Eqs. (4)-(8). Furter, since different 
SSCs, which may be regarded as different ways to specify 



center-of-motion world line of each spinning body, can be 
related by a PN order shift of the center-of-motion world 
lines [see Appendix A of Ref. [7]], it is clear that the 
equations of motion, given in Ref. [7] and that given by 
Eq. (1), for the two cases considered, describe the same 
dynamics. More importantly, using Eqs. (2.6)-(2.13) of 
Ref. [7], it is straightforward to show that the quantities 
E, L, S, and L ■ S are also conserved to the required PN 
order for our two cases. This shows that the determinis- 
tic nature of the dynamics of spinning compact binaries, 
for the two distinct cases (i) and (ii), does not depend 
on coordinate conditions, SSCs, and the representations. 
It is interesting to note (again) that Ref. [11] reported 
chaos in case (i), where only one object spins, using the 
dynamics prescribed in Ref. [7]. 

The physical argument employed in Ref. [11] to explain 
the apparent "chaos" is the highly irregular behavior of 
the orbital plane due to the inclusion of leading order 
spin effects on to the conservative 2PN accurate dynam- 
ics, associated with the nonspinning compact binaries. 
In Ref. [11], the Newtonian orbital angular momentum 
vector was used to specify the orientation of the orbital 
plane. We believe that it is not an appropriate quantity 
to specify the orbital plane as it is not even perpendic- 
ular to the orbital plane, when spin effects are included. 
However, in the Hamiltonian formulation of PN dynam- 
ics, L = r x p, always remains perpendicular to the or- 
bital plane and easily tracks its orientation. Recall that 
L, as defined above, also occurs in the Poincare algebra 
of gravitating two-body system of spin-less particles to 
3PN order [23]. We observe that in the two cases, where 
analytic solutions exist, the orbital plane, specified by L, 
carves out a simple cone as it precesses around the invari- 
ant direction, defined by the total angular momentum J. 
Further note that, in these cases, L = —S and we think 
this is not the way the evolution of L and S is treated 
in Ref. [11] [see Eq. (2.6) in Ref. [11]]. When spin-spin 
effects are included, the angle between L and J no longer 
remains a constant, as is roughly evident from Eq. (4.15) 
in Ref. [13]. However, we note that the time scale for the 
evolution of 9, which defines the angle between L and J, 
should be large and the amplitude of its variation small. 
This is so as spin-spin interactions, responsible for the 
variation of 9, may be treated to be 2PN or 3PN order 
corrections, depending on the magnitudes of spins [see 
discussions at the end of Sec. II in Ref. [13]]. 

It was stated in Ref. [11] that beyond 2PN order, the 
dynamics of compact binaries will be chaotic even in the 
absence of spin. The erroneousness of this conclusion 
is clearly evident from the recent determination of 3PN 
accurate generalized quasi-Keplerian parametrization for 
the solution of the 3PN accurate equations of motion for 
two nonspinning compact objects moving in eccentric or- 
bits [12]. The underlying dynamics, derivable from an 
ordinary Hamiltonian, can be extracted from the papers 
of Damour, Jaranowski, and Schafer [24]. It can easily 
be shown that the motion is restricted to a plane, namely 
the orbital plane, and we can introduce polar coordinates 
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such that r = r(cos tp, sin tp, 0). The 3PN accurate gen- 
eralized quasi-Keplerian parametrization for compact bi- 
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I = n (t — to) 
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where v 



The PN accurate orbital elements (Xf , ^/-i ^5 6^ i ^5 Enid 
and the PN orbital functions g 4t , g 6t , / 4t , fet, i&t, Kt, 
hip, fs v , 9iip, 96ip, H v , and h 6ip , expressible in terms of 
E, L, and 7/, are obtainable from Ref. [12]. The number 
of degrees of freedom are three and there are sufficient 
conserved quantities to restrict the orbits to a plane and 
the dynamics to be integrable. The deterministic na- 
ture of the dynamics does not change when radiation 
reaction, which enters at 2.5PN order, is included. It 
is actually possible to incorporate, almost analytically, 
the effects of reactive -ff2.5PN(»",p), derivable easily from 
Ref. [14], in to the 3PN accurate conservative dynam- 
ics of Ref. [12], employing the techniques described in 
Ref. [22]. We would like to remind that a fully 2.5PN ac- 
curate orbital dynamics of inspiralling binaries consisting 
of nonspinning compact objects of arbitrary mass ratio, 
available in Ref. [22], requires numerical evaluation of 
just one quadrature [see Sec. V of Ref. [22] and discus- 
sions therein]. Further, with the help of Refs. [12, 22], 
one can easily conclude that the motion is still restricted 
to a plane to 3PN order in the case of inspiralling non- 
spinning compact binaries. 

Finally, note that PN accurate differential equations, 
employed in Ref. [11], are realistic, only if one is allowed 
to neglect the omitted higher order correction terms. 
Therefore, extreme care should be taken to make sure 
that the solutions to these differential equations are also 
as PN accurate and this is not an easy task, while solving 
these differential equations numerically. In the paramet- 
ric solutions, described above, we can easily control their 
PN accuracy and determine the validity of the approxi- 
mate dynamcis in a gauge invariant way. However, it is 
not always possible to obtain PN accurate simple para- 
metric solutions, as given in Refs. [12, 13]. For example, 
in case of a spinning compact binary with leading order 
spin-orbit interaction, it will not be possible to obtain 
a simple parametric solution for the general case, where 
mi 7^ in 2 and Si ^ S 2 , as is no longer a constant an- 
gle there. In this case, it is still possible to parametrize 
the radial motion, but the angular motions do not ac- 



naries moving in eccentric orbits, in the center-of-mass 
frame, is then given by 

(9a) 
(9b) 

(9c) 
(9d) 

I 

commodate any simple parametric solutions. 

We want to point out that our case (ii) , which includes 
two spinning objects with leading order spin-orbit inter- 
action, should not be taken to argue that spin-spin in- 
teractions will not make the dynamics of spinning com- 
pact binaries chaotic. Our case (ii) should be viewed as 
a toy model to test numerical investigations that probe 
the deterministic nature of PN accurate dynamics of spin- 
ning compact binaries. When Hss terms, which specify 
the spin-spin interactions, are included in Eq. (1), even 
the radial motion does not allow simple parametrization. 
However, it should be possible to obtain semi-analytic 
PN accurate solution to the PN accurate dynamics of 
spinning compact binaries, moving in eccentric orbits, 
when all leading order spin effects are included. This will 
require adapting techniques from classical perturbation 
theory and will be similar to the way classical spin-orbit 
coupling was incorporated into the Newtonian accurate 
orbital motion [25, 26]. We believe that such PN accu- 
rate solutions to PN accurate binary dynamics may be 
more suitable to analyze the question of "chaos" in PN 
accurate dynamics of spinning compact binaries, in a re- 
alistic manner. Many of these issues are currently under 
investigation in Jena. 



III. SUMMARY 

Let us recapitulate. We have shown that the conserva- 
tive PN accurate dynamics of spinning compact binaries 
with leading order spin-orbit interaction, for the two dis- 
tinct cases mentioned above, is integrable and hence can- 
not be chaotic. We emphasized that these two scenarios 
are as deterministic as the PN accurate conservative dy- 
namics of nonspinning compact binaries. We hope that 
our parametric solutions, for the two distinct cases, will 
be employed to test numerical investigations that probe 
the question of chaos in spinning compact binaries. Fi- 
nally, we feel that PN accurate semi-analytic solutions 
to the PN accurate dynamics of spinning compact bina- 
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ries may help to clarify the question of chaos in spinning 
compact binaries. 
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